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Abstract 

A semiclassical, four-level model of a nanosecond, chirped pulse interacting with all optically 
accessible hyperfine states in the ultracold Rb atom is analyzed aiming at population inversion 
within 55 i/ 2 electronic state. The nature of two-photon adiabatic passage performed by such a 
single pulse having a bandwidth smaller than the hyperfine splitting of 5S'i/2 state is investigated 
in the framework of the dressed state picture. It is shown that two dressed states are involved in 
the adiabatic dynamics of population inversion. The excited state manifold appeared to play an 
important mediating role in the mechanism of population transfer. 
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Adiabatic passage known to be a vital control technique for the manipulation of dy¬ 
namics in atoms and molecules and the preparation of predetermined superposition states 
in these systems. It is broadly implemented in control schemes at ultracold temperatures. 
Examples include the use of the pulse area solution for fast two-qubit phase gates [1], a 
creation of ultracold molecules via stimulated Raman adiabatic passage (STIRAP) [2], and 
the implementation of different combinations of chirped pulses to deterministically excite 
Rydberg states M. Here we discuss control of valence electron in ultracold Rb by a sin¬ 
gle nanosecond pulse whose bandwidth is much narrower than the transitional frequency 
between hyperhne levels of 5^S'i/2 state. The population inversion is achieved through the 
Raman transitions that involve hyperhne structure of 5^Pi/2 or 5^P3/2 states and with the 
aid of linear chirping of the pulse. The adiabatic solution is found through a systematic nu¬ 
merical analysis of the response of a four-level system, representing all optically attainable 
hyperhne states, on a broad variation of the held parameters. We have made a detailed anal¬ 
ysis of the dressed state picture to gain insight into adiabatic mechanism of the two-photon 
Raman transition by means of a single narrow-band nanosecond pulse, which reveals the 
involvement of two dressed states into adiabatic passage producing population inversion. 
These two dressed states form a subset owing to two energetically close hyperhne states 
of the transitional 5^Pi/2 or 5^P3/2 state. Because only one chirped pulse is implemented 
having the bandwidth narrower than the hyperhne splitting of 5^S'i/2 state, the excited state 
manifold plays the key role in the passage as a mediator, thus, distinguishing this approach 
from previous experiences. 

We consider the four-level system which takes into account all optically allowed transi¬ 
tions between hyperhne states belonging to 5^S'i/2 and 5^Pi/2 or 5^P3/2 states, Fig.Q. The 
Hamiltonian that describes the four-level system interaction with a single chirped nanosec¬ 
ond pulse reads 


Hint h 


A -f- ca43 -f- Oi(t — P) 0 — l‘2 — VLjiit) 

0 /\.-\-ijJi3-\-oj21-\-ot{t — T') —r2/j(t)/2 

-nR{t)/2 0 0 

—flR(t)/2 —flR(t)/2 0 Ci;43 


( 1 ) 


Here flR{t) = —fiEQ{t)/h is the Rabi frequency with the peak value Qr, a is the linear 
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chirp parameter, a/27r has units Hz/s, and A is the one-photon detuning. Solving numeri¬ 
cally the time-dependent Schrodinger equation with the Hamiltonian in Eq.Q for various 
values of the held parameters provides an accurate picture of light-matter interactions and 
allows for hnding the exact values of the held parameters required to obtain a predeter¬ 
mined non-equilibrium superposition state, population invention or population return. It 
also reveals the adiabatic region of population transfer to the target state |2 >, which is 
the upper hyperhne state F=2 of the 5^S'i/2. Populations of the four states at the end 
of the pulse as a function of the pulse chirp rate and the full width at the half maximum 
(FWHM) are presented in Fig.(|^, [B]. FWHM of the Gaussian pulse relates to the pulse du¬ 
ration To as FWHM=ro2A/ln 2, ojq. The adiabatic region of light-matter interaction leading 
to population inversion is observed for parameters that satisfy the adiabaticity conditions 
|a/(27r)|ro > a; 2 i and |Q!/(27r)| < The physical values applicable to ®®Rb are, e.g., the 
peak Rabi frequency Qr = UJ 21 (3.035 GHz), the chirp rate a/27r=-0.3[ci;2i] (-3 GHz/ns) 
or faster and the pulse duration tq > 5.5a;^/, (> 1.8 ns). The negative value of the chirp 
rate is well understood since we start from the one-photon blue detuning with the largest 
transition frequency UJ 41 and gradually decrease the instantaneous frequency to pass through 
each one-photon resonance, hrst with LJ 41 , then with cusi, LJ 42 and 0732 . The exemplihed held 
parameters may be obtained in modern experimental setups such as described in, e.g., [7]. 
Since the spectral bandwidth of the nanosecond pulse ( 0.5 GHz for 1.8 ns pulse) is much 
narrower than the energy separation between the hyperhne states (3.035 GHz) of the 5^S'i/2, 
a question of fundamental interest arises as to what is the mechanism of the adiabatic popu¬ 
lation transfer performed with the two photons that are never present in the system with the 
frequency ’’right” to satisfy the two-photon resonant condition? We performed the dressed 
state analysis to gain insight into the adiabatic and nonadiabatic nature of quantum control 
of population dynamics in the four-level system using a single narrowband but chirped laser 
pulse. 

We hrst outline a basic concept of the dressed state analysis, pin], and its extension 
to the case when adiabatic passage may occur within a subset of dressed states coupled 
to each other. A wave function of a quantum system |T(f) > may be written as a linear 
superposition of the bare states in the held interaction representation \i > with the respective 
probability amplitudes Ci 
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|'I'(«) >= S*C.|i > , (2) 

The time-dependent Schrodinger equation then reads as ifiC = HintC. We apply a 
unitary transformation T to the Hint leading to diagonalization of the Hamiltonian. Here, 
T is an eigenvector matrix of Hint- The obtained Hamiltonian is the so called dressed 
state Hamiltonian Hd = H'^HintH written in the basis of the dressed states |/ > such that 
|'l/(t) >= > and Cd = TC. Then, putting the reverse expression C = T+C^ into 

the Schrodinger equation and assuming that all quantities are time dependent we arrive at 

ihCd = HdCd - Cd- (3) 

Since the Hamiltonian Hd is diagonal, the dressed states would evolve without mixing with 
each other if to disregard the second term on the right side. This is the essence of the 
adiabatic approximation, when the system, once placed in a selected dressed state by the 
initial conditions, continues evolution within this dressed state only. The second term is 
responsible for the nonadiabatic coupling between the dressed states, it contains matrix 
operator TT which is non-diagonal and determines the degree of non-adiabatic mixture 
between the dressed states. If the matrix elements of the TT are much less than the 
energy splittings between the respective dressed states, the dynamics may be considered 
as adiabatic. Analises of the time-dependence of the dressed state energies and the wave 
functions as well as a comparison of non-adiabatic and adiabatic terms help to estimate 
the degree of adiabaticity and a possibility for quantum control. From another hand, if to 
aim to hnd the held parameters that provide the adiabatic solution, it is useful to move 
to the dressed state basis and within the adiabatic approximation hnd the held conditions 
and parameters for dynamics in a single dressed state. When implemented in the exact 
Schrodinger equation, these parameters may yield quasi-adiabatic behavior within the exact 
Schrodinger picture. 

In a multi-level case, the non-adiabatic coupling may be small for some dressed states, but 
signihcant for the others. Then coupled dressed states may create a subsystem within which 
the dynamics occurs adiabatically. For the four-level system described by the Hamiltonian 
in Eq.Q, two dressed states, |/ > and \III >, are coupled in the vicinity of the avoided 


4 


crossing and provide adiabatic passage, while two other dressed states stay intact: 
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Note, that a superposition of two active dressed states, |/ > and \III >, may be approx¬ 
imated by a single dressed state if to reduce the four-level system to a three-level A system 
by substituting two energetically close excited states by a single one. We will discuss this 
approximation in details below. In the framework of the three-level A system, the adiabatic 
dynamics takes place within a single dressed state having the energy dependence on time 
resembling that of a superposition state in the four-level system. 

We demonstrate the concept by analyzing the time dependence of the dressed state 
energies in the four-level system and the squares of the respective eigenvector elements 
Tij that show the time evolution of the population of the bare states within each dressed 
state. The parameters of the nanosecond chirped pulse that provide the adiabatic passage 
are chosen from the numerical solution of the Schrodinger equation, Fig.(|^. Here a broad 
adiabatic region leading to population inversion to the upper hyperhne state of the 5S'i/2 
state is observed for the parameters starting from FWHM = 2.5 ns and higher and the 
absolute value of the chirp rate |Q;/27r| grater than 2 GHz/ns. 

The exact time-dependent picture of the adiabatic passage to the hnal hyperhne sate 
obtained by numerically solving the time-dependent Schrodinger equation with the Hamil¬ 
tonian in Eq.Q is shown in Fig.(|^ for parameters from the adiabatic region, e.g., FWHM= 
2.995 ns, a/ 27 i = -2.947 GHz/ns, 12^=3.035 GHz and A = 0. Population how from the 
ground |1 > to the excited |2 > state begins at about half-way before the Rabi frequency 
reaches the peak value. It is owing to the one-photon oh-resonance contribution of light- 
matter interaction into the population passage. Notably, population inversion dynamics 
follows the excitation of states |3 > and |4 >; those states get transitionally populated to up 
to 10% and show oscillations. The oscillations are attributed to small nonadiabatic coupling 
between all dressed states. They vanish with the increase of the Rabi frequency. 

The respective dressed state energies for the same held parameters FWHM= 2.995 ns, 
a/ 27 i = -2.947 GHz/ns, r2ii;=3.035 GHz and A = 0 are depicted in Fig.Q. It demonstrates 
that if the chirp rate is large enough, the change of the instantaneous carrier frequency is 
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FIG. 1: Four optically attainable hyperfine states of 5S and 5P shells with the energy differences 
that correspond to the D1 line. Initially, the population is in the ground state |1). Note that the 
hyperfine splitting of the 5Si/2 orbital is approximately an order of magnitude greater than the 
splitting of 5Pi/2 orbital, [TU] . 

sufficient to swipe adiabatically through the two-photon resonance. The dynamics in the 
system begins within the blue dressed state \I >, which coincides with initially populated 
bare state |1 >, (shown in blue dashed color). Further on, the blue dressed state |/ > 
approaches the green dressed state \III > near the peak value of the held amplitude to 
form an avoided crossing. In its vicinity, population moves efficiently from the blue |/ > 
to the green \III > dressed state and, thus, resides transitionally on the excited bare state 
manifold for a restricted period of time before moving to the hnal bare state |2 > as the 
green dressed state \III > evolves to become 100% constituted of it. Thus, the excited 
states 13 > and |4 > keep population for the time needed for the instantaneous frequency 
to acquire the value needed to accomplish the two-photon resonance. In such a way, the 
adiabatic passage is performed by a subset of coupled through the avoided crossing dressed 
states and requires the excited state manifold to mediate the dynamics owing to only one 
chirped, narrowband pulse used to perform the inversion. 

For a comparison, the case when the chirp is not large enough to provide the range of 
frequencies needed to satisfy the two-photon transition is shown in Fig(|^. Here the chirp 
rate is a/2Tr = -0.092 GHz/ns and the pulse duration is tq = 1.799 ns (FWHM = 2.995 
ns), giving aro= 0.166 GHz <C UJ 21 = 3.035 GHz. The other parameters are the peak Rabi 
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FIG. 2: The end-of-pulse population distribution in the four-level system, achieved via two-photon 
transitions using a single, linearly chirped laser pulse. The values of the system parameters are uj 2 i 
= 3.035 GHz, oji'i = 0.362 GHz, characteristic for [lO], the peak Rabi frequency is 12^=3.035 
GHz, and one-photon detuning A is zero. 

frequency VLr = a; 2 i=3.035 GHz, and one-photon detuning A = 0. Here, blue dressed state 
initially coincides with blue bare state |1 >; this picture remains for most of the pulse 
duration. Then, close to the exponential end of the pulse amplitude, blue dressed state 
approaches the green dressed state, which is mainly a superposition of the excited states 
13 > and |4 >, and nonadiabatically transfers a fraction of population to them before the 
pulse ceases. As the result, the population remains mostly in the initial, ground state owing 
to the held parameters not satisfying the adiabaticity condition for population inversion 
determined by the condition a/(27r)ro > UJ 21 , even though the Landau-Zener condition is 
satished, VL\/\a\ ~ 3/2 x 10^. 

To beneht from an analytical solution, we performed a comparative analysis with an 
effective three-level A system, Fig.([^ that works as a good approximation to the four-level 
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FIG. 3; Time-dependent picture of the population dynamics in four hyperfine states demonstrating 
adiabatic passage to the final state |2 >. The field parameters are ^}r=3.035 GHz, FWHM=2.995 
ns and a/(27r)=-2.947 GHz/ns. 
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FIG. 4: Dressed state energies as a function of time that lead to adiabatic passage to the final state 
|2 >. The field parameters are Hi?=3.035 GHz, FWHM=2.995 ns and a/(27r)=-2.947 GHz/ns. 

system giving a qualitatively similar quantum yield when the pulse duration and the chirp 
rate satisfy |a/(27r) |ro uj^s, and VLr 3> 0143 . A detailed numerical solution for the three-level 
A system is discussed in mm 

From the Hamiltonian in Eq.Q, we may easily get the following set of equations for 
the probability amplitudes in the held interaction representation assuming the one-photon 
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FIG. 5: Dressed state energies as a function of time for the field parameters Djj = 3.035 GHz, 
FWHM = 2.995 ns and Q;/( 27 r) = -0.092 GHz/ns. No adiabatic passage occurs to the final bare 
state 12 > and the dynamics is preserved within dressed state \I >, which coincides with bare 
states 13 > and |4 > at the end of the pulse. 



FIG. 6 : A three-level A system considered as an approximate model with two optically attainable 
hyperfine states of 5P shell described by a single transitional state |3 >. The approximation is 
implemented owing to the hyperfine splitting of the f)Pi /2 orbital being about an order of magnitude 
less than the splitting of 55^/2 orbital. It is valid under the conditions that |a/(27r)|ro a .’43 and 

D/j S> W 43 . Initially, the population resides in the ground state |1). 
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detuning A=0 


icLi = a{t — T)ai — flR{t){a3 + 04)72 
id2 = [ui2i + a{t - T)]a2 - ilR{t){a3 + 04)72 

ids =—f 2 /j(t)(oi + 02)72 — Ci;4303 ( 5 ) 

id 4 = -(2i?(t)(oi + 02)72 


By making a substitution 


(03 + 04)7 a /2 = 0 + ( 6 ) 

(03 - 04 ) 7 a/ 2 = o_, (7) 

we arrive at the following relation for 0 + and o_ 

i(d+) = i(d 3 + 04 ) 7 a /2 = -f 2 R(f)(oi + 02 ) 7 a /2 - 60437( 8 ) 
i(d_) = i(d 3 - d4)7\/2 = -u^s/V^a^. (9) 


If U43-\/2 is small enough compared to it may be neglected. Then, 03 — 04 = const, 
and o_ may be omitted from the dynamics calculation. Finally, the Eqs. ([^ are reduced to 
a set of three coupled differential equations 


idi = a(t — T)oi — f2i^(t)7'\/2o+ (10) 

id2 = [6021 + ait - T)]o2 - flR(t)/^/2a+ (11) 

id+=-QR{t)/y/2{ai + 02 ) ( 12 ) 

The Hamiltonian for the three-level approximation in the held interaction representation 
reads 


Hint h 


a{t-T) 0 -VtR{t)/\/2 

0 6 O 21 -|-cr(t — T) —VLR{t)/\/2 

-nR{t)/V2 -^R{t)/V2 0 


(13) 


Analytical diagonalization of the Hamiltonian in Eq.(13) leading to expressions for the 
dressed state energies and respective eigenfunctions showed no dark state solution as it was 
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the case, e.g., in conventional STIRAP scheme, [13]. The time-dependent wave function 
describing each dressed state contains nonzero probability amplitudes for all three bare 
states. Since the expressions for the dressed state energies and the probability amplitudes 
look heavy, we do not present them here, but rather show their time dependence obtained 
numerically. A numerical analysis of the dressed states was performed for the three-level 


A system described by the Hamiltonian in Eq. (13) within the same range of parameters as 
for the four-level system. As an example, the results for a/27r = -2.947 GHz/ns, FWHM 
= 2.995 ns and Qr = 3.035 GHz are discussed in more details. The time-dependence of 
the dressed state energies is depicted in Fig.Q. Here the dynamics occurs within a single 
dressed state |J >, shown in blue color. Initially, it coincides with bare state |1 > (dashed 
blue), followed by adiabatic transition within the same dressed state in the vicinity of the 
peak values of the pulse amplitude from the bare state |1 > through the excited state |3 > 
(dashed green) to hnal bare state |2 > (dashed black). In this approximate three-level model, 
the dressed state |/ > plays the role of the subset of two dressed states, |/ > and \III >, 
in the four-level case. In Fig.(|8|), the adiabatic dynamics of population transfer between the 


bare states within each dressed state is shown by the time dependence of T^y It reveals a 
smooth adiabatic passage from the initial bare state |1 > to the hnal bare state |2 > in the 
dressed state |/ >,Fig.(|^). 

However, if to choose a set of parameters such that they do not satisfy the condition 
Q;/(27r)ro > 0021-1 which is, for example, a/27r = -0.092 GHz/ns, FWHM = 2.995 ns, the 
dynamics still occurs within a single dressed state but does not lead to the population 
inversion to the hnal bare state |2 >. In the beginning, the energy of dressed state |/ > 
coincides with the bare state |1 >; however, the pulse ceases before the dynamics within 
dressed state |/ > progresses to the hnal bare state |2 >, which makes it retaining in the 
intermediate state |3 >, Fig.([^. The chirp rate of the pulse is not fast enough to switch 
through the two-photon resonance. The bare state population dynamics within each of 


three dressed states is shown in Fig. (10), it supports this outcome by manifesting adiabatic 
passage from the initial bare state |1 > to the intermediate bare state |3 > within the active 


dressed state |/ >, Fig.(lOr). 

In summary, we have analyzed the mechanism of population dynamics in a new method for 
population inversion within the hyperhne structure in alkali atoms at ultracold temperatures 
that relies on implementing a single nanosecond chirped pulse of intensity on the order of 
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FIG. 7: The adiabatic dynamics of population transfer between the bare states within each dressed 
state in the three-level A system. The field parameters are nij=3.035 GHz, FWHM = 2.995 ns 
and Q:/( 27 r) = -2.947 GHz/ns. 

kW/cm^ having the bandwidth much narrower than the two-photon transition frequency in 
the atomic system. The results are based on the developed semiclassical model of the pulse 
interaction with the four-level system representing all optically accessible hyperfine states 
of 5 ^S'i/2 and b‘^Pi /2 or b‘^P ^/2 states in ultracold ®^Rb. The adiabatic passage leading to 
population inversion is achieved for parameters that satisfy the condition |a/(27r)|ro > a; 2 i 
and the Landau-Zener adiabaticity condition |a/(27r)| < Dressed state analysis was 
performed to gain understanding about the mechanisms of two-photon Raman transitions 
performed by a single, narrowband, chirped pulse having the bandwidth Aw W 21 . It 
revealed an existence of a subset of dressed states coupled in the vicinity of avoided crossings 
that perform the adiabatic passage leading to the population inversion. When considered as 
an approximation of the four-level system, the three-level A system demonstrates population 
inversion within a single dressed state. This dressed state resembles the time-dependence of 
two dressed states in the active subset in the four-level system. This justifies the validity of 
several dressed states approach to adiabatic passage. 

The authors acknowledge fruitful discussions with Vladimir Malinovsky, Elena 
Kuznetsova and Phillip Gould. This research was supported by the National Science Foun- 
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FIG. 8: The adiabatic dynamics of population transfer between the bare states within each dressed 
state. Population inversion between states |1 > and |2 > is observed within the dressed state |/ >, 
(a). The field parameters are = 3.035 GHz, FWHM = 2.995 ns and Q;/(27r) = -2.947 GHz/ns. 
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FIG. 9: Dressed state analysis of energy picture for the three-level. The field parameters are 
= 3.035 GHz, FWHM = 2.995 ns and a/(27r) = -0.092 GHz/ns. 
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FIG. 10: The adiabatic dynamics of population transfer between the bare states within each 
dressed state which does not lead to population inversion between states |1 > and |2 >. The field 
parameters are Qr = 3.035 GHz, FWHM = 2.995 ns and a/(27r) = -0.092 GHz/ns. 
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